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ABSTRACT 


The  intensity  of  second  harmonic  light  generated  by 
reflection  frcia  a  Gall inn  Arsenide  crystal  is  calculated  in 
detail  in  four  cases  of  polarization.  Relevant  properties  of 
GaAs  are  stated  and  an  experiment  is  described.  Based  on  this, 
constitutive  relations  are  written  in  the  weak  field  approxi¬ 
mation  and  the  second  and  third  order  material  tensors  are 
reduced  using  crystal  symmetry.  The  nonlinear  wave  equation 
is  solved  in  the  parametric  approximation  when  phase  matching 
is  not  present.  The  boundary  value  problems  for  the  funda¬ 
mental  entering  the  medium  and  for  the  second  harmonic  leaving 
are  solved.  The  only  significant  null  solution  is  found  to 
be  the  Brewster’ s  angle  extinction.  Seme  numerical  exploration 
of  the  exact  effects  of  absorption  at  both  fundamental  and 
second  harmonic  are  made.  They  are  found  to  be  small  and  are 
negligible  for  qualitative  purposes  unless  the  imaginary  part 
of  the  dielectric  constant  exceeds  the  real  part.  Variation 
of  intensity  with  crystal  orientation  is  described  in  detail. 
The  occurrence  of  phase  matching  as  a  null  solution  of  the  wave 
equation  is  seen  to  offer  a  natural  approach  to  the  problem 
of  obtaining  the  geometry  of  phase  matching  for  the  most 


general  optical  medium. 
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SECOND  HARMONIC  GENERATION 
BY  INFLECTION  FROM  GaAa 


1.  THE  GENERAL  PHASE  MATCHING  PROBLEM 

The  onset  of  stimulated  effects  in  nonlinear  optics  is  quite 
dramatic.  Ordinary  nonlinear  effects  arising  from  even  the  most 
intense  excitation  are  barely  detectable  whereas  stimulated  effects 
may  be  of  the  same  order  of  intensity  as  the  excitation.  It  may 
even  be  said  that  the  extension  of  stimulated  effects  into  the 
optical  region  is  the  cause  of  the  recent  explosion  of  interest 
in  nonlinear  optics.  For  the  nonlinear  effect,  second  harmonic 
generation,  a  criterion  of  analogous  importance  is  that  of  phase 
match  between  the  fundamental  and  second  harmonic.  When  this  con¬ 
dition  is  met,  the  two  waves  proceed  in  the  same  direction,  in  step, 
in  principle  allowing  complete  power  transfer  from  the  fundamental 
beam  to  the  second  harmonic  beam.  If  thi6  condition  is  not  met, 
the  second  harmonic  intensity  remains  very  small  even  if  a  very 
intense  fundamental  is  used. 

Often  a  medium  is  normally  dispersive  and  the  index  of  re¬ 
fraction  at  the  second  harmonic  is  higher  than  at  the  fundamental. 
Then  it  is  not  possible  for  the  two  waves  to  proceed  in  the  same 
direction  and  be  in  6tep  for  the  magnitude  of  their  wave  vectors 
is  different.  The  most  well  known  way  of  achieving  phase  match  is 
to  use  the  birefringence  of  certain  anisotropic,  nonabsorbing, 
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nonlinear  crystals.1  In  an  anisotropic  medium  there  are  in  general 
two  indices  of  refraction  for  each  direction  and  they  vary  with 
direction.  Thus  sometimes  a  direction  can  be  found  where  an  index 
at  the  fundamental  is  the  same  as  at  the  second  harmonic.  Recently 
the  use  of  optical  activity  has  been  advanced  as  a  means  of  achieving 
phase  match.2  Both  of  these  methods  make  use  of  linear  optical 
properties,  which  may  be  comprehensively  categorized  by  the  use  of 
matrix  algebra.  If  A’  is  the  2x2  matrix  obtained  by  projecting 
A  in  E  =  AD  onto  the  plane  transverse  to  the  wave  vector,  then 
the  1)  real  symmetric,  2)  real  antisymmetric,  3)  imaginary  symmetric, 
4)  imaginary  antisymmetric  parts  of  A’  contribute  respectively  the 
following  optical  effects:  l)  birefringence,  2)  optical  activity, 

3)  anisotropic  absorbtion,  4)  circular  dichoism.3  In  general  optical 
media  exhibit  all  combinations  of  these  effects.  From  this  it  is 
clear  that  if  the  two  methods  above  represent  fairly  the  state  of 
results  in  the  general  phase  matching  problem,  as  the  author  believes, 
then  it  is  not  completely  solved. 

This  report  is  a  calculation  of  the  second  harmonic  intensity 
produced  by  reflection  from  a  Gallium  Arsenide  crystal,  including 
the  effects  of  absorbtion  at  both  fundamental  and  second  harmonic. 

1  D.  A.  Kleinman,  Phys.  Rev.  128,  1761  (1962). 

2  H.  Rabin  and  P.  Bey,  Bull.  Am.  Phys.  Soc.  12,  8l  '1967). 

3  G.  Ramachandran  and  S.  Ramaseshan,  Handbuch  der  Physik  XXV/l 
107  (1961). 
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Conceived  of  as  an  approach  to  the  general  phase  matching  problem, 
second  harmonic  production  by  reflection  is  rather  oblique.  In 
fact  it  is  not  even  a  true  subcase  since  the  parametric  approxi¬ 
mation  used  here  is  not  capable  of  eliciting  the  effects  at 
phase  match.  Still  it  exhibits  many  necessary  techniques  and  a 
thorough  understanding  of  it  is  no  doubt  a  prerequisite  for  an 
attack  on  the  general  phase  matching  problem.  Thus  while  most 
of  the  theory  appearing  here  has  appeared  before,4  general  methods 
are  used  where  known  because  of  the  above  outlook  on  the  general- 
phase  matching  problem. 

The  essential  steps  in  the  calculation  of  the  reflected 
second  harmonic  intensity,  each  giving  rise  to  a  section  below, 
are:  describing  GaAs  and  the  geometry,  solving  the  nonlinear 
wave  equation,  solving  the  boundary  value  problems  and  obtaining 
the  intensities.  The  next  to  last  section  gives  some  numerical 
results  and  the  last  section  remarks  on  the  general  phase  matching 
problem. 


4  P.  Butcher,  Nonlinear  Optical  Phenomena,  Ohio  State  University 

(1965) 
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2.  GALLIUM  ARSENIDE  AND  THE  EXPERIMENTAL  SETUP 
An  experiment  with  GaAs  of  the  sort  considered  here  was  first 
reported  in  196l.5  Since  GaAs  is  opaque,  a  reflection  experiment 
appears  to  be  required.  It  has  point  group  symmetry.6 7  This 
point  group  is  cubic  so  second  order  material  tensors  reduce  to 
scalars  and  linear  optical  effects  become  isotropic  7  Furthermore 
the  point  group  excludes  the  antisymmetric  part  of  a  second  order 
tensor  so  optical  activity  and  circular  dichroism  are  absent.  Thus 
the  linear  optical  properties  of  GaAs  may  be  completely  described 
by  a  complex  dielectric  susceptibility 

e  =  (v  +  1  h)2  (2.1) 

where  v  and  h  are  the  conventional  index  of  refraction  and  absorbtion 
coefficient.  Suppose  the  illumination  is  ruby  laser  light  at  6943A. 

If  the  corresponding  radian  frequency  is  uoQ,  then  for  present  purposes 
the  linear  optical  properties  are  completely  described  by8 

5  J.  Ducuing  and  N.  Bloembergen,  Phys.  Rev.  Letters  10,  '+74  (1961). 

6  Wycoff  Structure  Tables  Vol.  I 

7  The  latter  part  of  this  section  is  a  discussion  of  the  effects  of 
symmetry  on  material  tensors. 

8  J.  Davey  and  T.  Pankey,  J.  Appl.  Phys.  35,  2203  (1964).  The  accuracy 
of  these  figures  is  probably  a  few  percent  of  the  largest  figure. 

They  were  obtained  by  a  Kramers-Kronig  analysis  of  reflectance  d-:.ca 
of  H.  R.  Philipp  and  H.  Ehrenreich,  Phys.  Rev.  129,  1550  (1963)- 
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Fig.  3  -  Crystal  Orientation 
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.25 


v  =5.79 

1 

v  =3-89 

2 

where  v  is  written  for  v(rxu  ). 
n  o' 

Specifically  the  experiment  considered  is  the  following. 

A  GaAs  crystal  is  prepared  with  a  (110)  type  plane  as  a  reflecting 
surface  in  such  a  way  as  not  to  damage  the  crystalline  properties 
of  the  surface  layers.  This  is  important  since  the  absorption 
depths  are 


h  — 
l 


k  =  1.53 
2 


(2.2) 


X/h  2:  4400A  ,  X/k  =  36lA. 

1  2 

Ruby  laser  light  of  known  intensity  I  and  polarization,  normal 
n  or  parallel  p,  strikes  the  reflecting  surface  with  angle  of 
incidence  0.  The  reflected  beam  is  filtered  to  remove  the  funda¬ 
mental.  The  second  harmonic  intensity  is  then  measured  in  two 
polarization  components  Jn  and  Jp.  The  crystal  may  be  rotated 
about  a  perpendicular  to  the  surface  to  expose  different  crystal 
orientations.  This  orientation  angle  is  Y  and  the  ¥  =  0  reference 
will  be  specified  later.  These  points  are  illustrated  in  Figs.  1-3* 
From  the  foregoing  and  from  general  considerations  one  may 
write  constitutive  relations  taylored  to  GaAs  and  the  experiment. 
For  a  medium  at  optical  frequencies 
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E*(ui)  =  E(-tu). 


(2.6) 


For  E(  t)  and  D( t)  simultaneously  real 

X^*{u))  =  x(l)(-)>  (2*7) 

X^*(u)  ,  u)  )  =  x^(-w .  j  -«  )•  (  2-8) 

12  1  2 

The  discussion  will  now  be  restricted  to  the  steady  state 
case  where  frequencies  present  at  one  time  are  always  present. 

Let  the  incident  beam  be  of  the  form 

E(  t)  =  Ee  lU)ot  +  cc, 

where  cc  represents  the  complex  conjugate  of  the  preceeding 
term.  In  the  frequency  domain, 

E((ju)  =  2nj(uu  -  u)q)E  +  2nj(u)  +  ujq)E* 

=  2tt5(u)  -  u)q)E  +  cc  (2.9) 

where  cc  implies  that  (2.9)  satisfies  (2.6).  Substituting  in 
(2.5) 
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D(u))  =  2tt6(uj  -  u)Q)[l  +  )]E  +  cc 

+  2rr6(u)  -  2u)q)x/ 2^(u3q,  “>0)E2  +  cc 
+  2tt6(u))x^(<u0>  -ouq)EE*  +  cc, 

using  (2.7)  and  (2.8).  Subsequently  the  u>  =  0  component  will 
be  ignored.  Reverting  to  the  time  domain,  holding  tuQ  constant, 

D(t)  =  e(u>o)Ee"lu,ot  +  cc  +  IhtP^  auQ,  t)  (2.10) 

where  e(u>o)  =  1  +  X^^)  (2.11) 

is  the  usual  susceptibility  and 

W^o’  *  X1Jk(2,(V  “c,)EAe’2‘'”0t  +CC  (2-12) 

Summation  on  repeated  indices  is  implied.  Henceforth  NL  will 
be  omitted  since  the  linear  polarization  is  not  mentioned.  Evi¬ 
dently  the  incident  wave  excites  second  harmonic  polarization 
of  the  medium  and,  as  will  be  seen  in  the  next  section,  second 
harmonic  radiation  is  emitted. 


9 


It  has  been  mentioneu  that  c  reduces  to  a  scalar  because  of 

(2) 

crystal  symmetry.  Likewise  the  second  order  susceptibility  x  > 
having  at  most  27  independent  components,  undergoes  considerable 
reduction.9  As  an  example  of  reduction  due  to  symmetry,  consider 
the  syumetry 


*ljk2)(w0-  “o'  ■  xlkJ<2)(‘"o-  “o'* 


( 2) 

obtained  from  (2.12).  The  9  elements  of  the  form  x..,  are 
unaffected  but  the  remaining  l8  are  equal  in  pairs.  Thus  it 
should  require  at  most  l8  constants  to  specify  X-y^  ' •  Applying 
the  elements  of  to  the  remaining  elements  reduces  them  to 
6  equal  elements  while  the  remainder  vanish.  There  are  several 
listings  of  tensors  reduced  by  crystal  symmetry  in  the  literature. 
One  such  for  third  order  tensors  gives  S  =  S  =  S  .10  In  this 

14  S  36 

"contracted"  notation  one  writes  the  components  symmetric  in 
the  last  two  indices  separately  from  those  antisymmetric.  Under 
this  convention  the  second  index  means  l:xx,  2:yy,  5:zz,  4:yz, 
5:xz,  6:xy.  Thus  in  present  notation 


133 


XP(  1®) 


9  A  general  treatment  of  this  problem  is  P.  Erdos,  Helv.  Phys. 
Acta  27,  495  (1964). 

10  J.  Giordmaine,  Phys.  Rev.  1^8,  Al601  (1965). 
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where  P(-125)  is  any  permutation  of  125  •  Henceforth  x  represents 

*133 .  This  redaction  is  performed  in  a  special  coordinate  system 
in  which  a  simple  representation  of  the  symmetry  elements  is 
possible.  For  many  point  groups,  for  cubic  point  groups  in 
particular,  the  axes  of  this  special  system  are  the  same  as 
those  of  the  conventional  unit  cell.11  The  result  is  for 
(2.12),  writing  the  complex  amplitudes  of  the  positive  frequency 
parts  only 


P 

l 


XE  E 

2  3 


> 


P  = 
2 


XE  E 

*  1  3 


> 


P 

3 


XE  E  , 
*  l  s 


a  local  relation  in  crystal  coordinates. 

This  relation  will  now  be  evaluated  in  experimental  coordi¬ 
nates  by  substituting  components  arising  in  two  rotations.  The 
rotation  equations  are  written  as  substitutions,  i.e., with  old 
coordinates  on  the  left, 


x  =  Rx'  . 


11  International  Tables  for  X-Ray  Crystallography,  Kynoch  Press, 
Birmingham,  England  (1952)  gives  the  convention  for  associating 
coordinate  axes  with  conventional  unit  cells. 
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In  this  way  one  may  drop  the  primes  immediately  as  there  is  nr, 
further  reference  to  the  old  coordinates.  Let  the  first  rotation 
he 


since  then  the  new  Z  axis  is  a  [110]  type  axis.  This  gives 


P  =  -  x[E  E  (cos2Y  -  sin2Y)  +  (E  2  -  E  2)sinYcosY]  cosY 
1  1  2  1  2 

+  |  v(E  2  -  E  2cos2Y  +  2E  E  sinYcosY  ••  E  2sin2Y)  sinY 

3  1  1  2  2 
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P  =  x[E  E  (cos2Y  -  sin2*)  +  (E  2  -  E  2)  sinYcosY]  sinY 
2  1  2  1  2 

+  i  x(®  2  *  E  ^os2*  +  2E  E  sinYcosY  -  E  2sin2Y)cosY 

31  12  2 

P  =  yE  (E  sinY  +  E  cosY) . 

3  3  1  2 

The  plane  of  incidence  is  the  XZ  plane  and  Y  is  the  angle  from 
the  ^ 100 )  plane  to  the  plane  of  incidence.  It  is  convenient 
to  reduce  this  to  cases  of  incident  field  polarization  normal 
and  parallel  to  the  plane  of  incidence. 


P 

1 
P 

2 

P 

3 


\  2sinY(3cos2Y  -  1) 
2 

-  3  XE  2sin2YcosY 
^  2 

0 


E  (E  =  0  =  E  ) 
n  i  3 


(2.13) 


P 

1 

P 

2 
P 

3 


-  |  xEi2sinYcos2Y  +  |  xE32sinY  J 

\  xEi2cosY( 3  sin2Y  -  1)  +  £  xE  2c°sY>  E^(E 
XE  E  sinY  1 

13  / 


0)  (2.14) 
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3.  THE  NONLINEAR  WAVE  EQUATION 

The  Maxvell  equations  for  media  at  optical  frequencies  are 

V  X  H  =  |  B,  7  •  B  =  0, 

(5-1) 

1  • 

7  X  E  =  -  -  D,  7*1*0. 

Applying  (2.3)  to  (3-1)  and  eliminating  B  and  H  gives  the  wave 
equation 


1  “ 

7  X  V  X  E  +  —  D  =  0, 

~  c2  ~ 

which  gives  in  the  frequency  domain 

iue 

7  X  7  X  E(<d)  -  —  D(»)  =  0.  (3.2) 

c2  ~ 

This  equation  is  nonlinear  since  D  is  nonlinear  in  E  so  the 
method  of  eigenfunctions,  in  which  the  general  solution  is 
obtained  by  superposition,  does  not  work.  Thus  E(u>)  must  be 
the  total  field.  One  anzatz  for  the  form  of  the  total  field 
is 

CO 

E(uj)  =  2tt^  6(uu  -  ou0)En  (3-5) 

n=-“> 
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1  ecause  once  one  frequency  is  introduced  the  sura  and  difference 
effects  of  (2.5)  may  build  up  all  harmonics.  Introducing  (3 .3) 
in  (3*2)  gives  an  infinite  system  of  connected  differential 


equations  for  which  the  incident  field  E  is  known.  An 
iteration  procedure  may  be  set  up  based  on 


where  r  is  the  order  of  the  iteration,  X  is  an  expansion  parameter 
to  be  set  equal  to  one,  and  a(n)  is  the  order  of  first  appearance 
of  the  n —  harmonic.  The  zeroth  iteration  is 


x  v  x  E  (°)  _  e  E  (0)  =  o 

1~1 


(3.*0 


where  eQ  =  e(nuo).  The  first  iteration  is 

2  2 

VX7X!^  .  k  -2-  e  E  W  =  l6rr—  P(E  (3.5) 

2  -2  2~2  „2 - 1 


and  there  will  be  an  u>  =  0  contribution  in  this  order  also. 

The  next  iteration  gives  a  correction  to  the  fundamental, 

( 2)  ( 2 ) 

E  ,  and  the  third  harmonic  appears,  E  '  .  A  correction 

~3 

to  the  second  harmonic  does  not  appear  until  the  third  iteration 

(2) 

and  it  is  third  order  in  x  •  a  very  small  quantity.  Thus 
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(3-5)  gives  an  accurate  E  .  This  approximation  is  called  the 

~2 

parametric  approximation  since  E  evidently  may  be  thought  of 

~i 

(2) 

as  a  sovrce  of  energy  which  generates  E  via  the  parameter  v  . 

~2 

If  significant  depletion  of  E  occurs,  the  approximation  is 

-1 

invalid.  Therefore  phase  matching  for  any  harmonic  actually 

present  must  be  excluded.  As  a  result  the  incident  wave  behaves 

as  it  would  in  a  linear  medium  and  the  second  harmonic  is  obtained 

as  the  solution  of  a  linear  inhomogeneous  wave  equation  with  an 

effectively  known  source  term  P(E  ). 

~  —l 

Using  the  complex  amplitudes  of  the  positive  frequency 
parts  of  (2.12),  the  source  term  has  the  plane  wave  form 

eik(u>0)  •  x  E  e‘£Uo)  •  x  e‘k(u)0)  •  x 

i  *ijk  lJ  ik 

hence  k(<li)o)  =  2k(u)Q)  (3*6) 

0,14  pi  - 

Because  of  homogeneity  these  amplitudes  are  constant  throughout 
the  medium.  Variation  due  to  absorption  is  accomplished  by 
using  a  complex  wave  vector.  For  (3*5)  then12 

12  Note  that  vector  magnitudes  are  written  |k|  =  k. 
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2 
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(e  $  e  since  phase  matching  has  been  excluded.)  In  the  last 
1  2 

step  (3*6)  was  used  and  (3-10)  was  anticipated.  To  complete  the 
solution  any  solution  of  the  homogeneous  equation  may  be  added. 
Homogeneous  solutions  are 


k(ou)  •  E^(uu)  =  0, 


(3-9) 


(3.10) 
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4.  T!IE  BOUNDARY  VALUE  PROBLEMS 


In  the  coordinates  set  up  in  section  2,  let  the  boundary  be 


the  plane  n  •  x  =  0  with  fi  pointing  into  the  first  medium  which 
is  transparent  and  nondispersive.  The  second  median  16  described 
by  section  2.  Given  an  incident  uniform  plane  wave  with  linear 
polarization,  find  the  reflected  second  harmonic  intensity.  There 
are  three  steps.  Find:  l)  the  transmitted  fundamental,  2)  the 
second  harmonic  polarization  (section  2),  and  3)  the  reflected 
second  harmonic.  This  section  concerns  the  first  and  last  steps. 

Rays  are  designated  in  Table  1  and  Fig.  4.  The  first  and 
second  symbols  on  each  line  pertain  respectively  to  the  fundamental 
and  second  harmonic.  Angles  are  measured  frcm  the  vertical  in 
the  respective  median.  In  the  boundary  value  problem  for  the 
fundamental  first  the  wave  vectors  will  be  determined  before 
dealing  with  the  amplitudes. 


■1 

Ampli¬ 

tude 

Wave 

Vector 

Angle  of 
Incidence 

Incident  or 

E 

k 

e 

Source  Wave 

S 

_s 

c  , 

Transmitted 

F 

n 

<t> 

Wave 

ri 

H 

T 

Reflected 

G 

m 

7 

Wave 

R 

r 

p 

Table  1  Ray  Nomenclature 


From  (3.10) 


k»-» 

=  m  =  — 

c 


(4.1) 


=  Is  =  v  2  .  V' 2  +  2  i  V  •  l" 


(4.2) 


where  t  =  V  +  it" ,  V  and  V'  real.  Since  the  waves  must  be 


the  same  function  of  position  on  the  boundary, 

k  •  x  =  V  ‘  x  =  m  •  x,  n  •  x  =  0 


(4.3) 


0  =  V'  •  x,  n  •  x  =  0. 


The  latter  implies 


l"  =  -3 


(4.4) 


(4.5) 


V  •  C  =  Vl"*OG+ 

•%.  -w  w 


(4.6) 


since  the  transmitted  wave  is  nonuniform  due  to  absorption  and 
decreases  most  rapidly  in  the  direction  -ft.  Using  the  vector 


identity 


n  x(xxn)=x-n-xn 


x  =  ft  x  ( x  X  n) ,  ft  •  x  =  0, 


(4.3)  becomes 


kxn•xxn  =  't.,  xn*xxn  =  mxn*xxn. 


But  since  a  vector  of  the  form  A  X  n  is  in  the  surface  and  since 
x  xn  is  any  surface  vector. 


k  x  n  -  xfi  =  mxn 


(4.7) 


ksin0  =  V  sin^»  =  msin7 


(4.8) 


Being  real,  m  is  now  completely  determined  by  the  magnitude 
(4.1),  the  two  components  (4.7)  and  by  knowing  in  which  medivm 
it  lies.  Then 


7  =  0. 


(4.9) 


Also  (4.7)  states  that  k,  V  ,  m,  and  n  lie  in  a  common  plane, 
defined  to  be  the  plane  of  incidence.  To  complete  the  deter¬ 
mination  of  £  define 


cJM 

n  i  - 
l  0) 

o 


j  «  =  cT 
i 


(4.10) 


Then  (4.2)  and  (4.6)  give 


e  =  n'2-n"2  +  2m  '  n 
ii  l  li 


and  using  (4.8) 


cos  =  - — V n  ' 2  -  sin20,  sine  <  n  '  , 
«  i  J-  i 


21 


taking  the  positive  square  root  since  by  (4.5),  V  and  lie 
in  the  same  quadrant.  This  gives 

c  =  n  ,2  -  n  "2  +  2*n  Vn  • 2  -  sin^8  (4.11) 

11111 

together  with  the  restriction 


sin0  <  n  ' 
l 


9 


(4.12) 


excluding  total  reflection.  The  inverse  of  (4.11)  is 


2n  ' 2  =  e  *  +  sin20  +  V(e  '  -  sin2©)2  +  e  "2 

ii  l  i 


2n  "2  =  -  e  1  +  sin20  +  V(e  1  -  sin20)2  +  e  1,2 
l  l  l  l 


(4.15) 


giving  the  "refractive  index"  and  "absorption  coefficient" 
in  terms  of  the  complex  susceptibility  and  angle  of  incidence. 
Since  the  definition  (4.10)  makes  n  '  and  n  "  dependent  on  0, 
they  are  not  strictly  material  constants  but  depend  on  a  non¬ 
material  boundary  condition  as  well.  The  usual  convention  is 
(2.1)  which  is  used  here  only  to  determine  e  from  (2.2).  Now  V 
and  V'  are  determined  by  (4.10),  (4.7),  and  (4.5)  just  as  m  was. 
Because  of  medium  homogeneity  this  determination  of  l  and  m 
applies  throughout  the  appropriate  medium. 
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The  amplitudes  may  be  determined  by  the  tangential  boundary 
conditions  for  electric  and  magnetic  vectors, 

(E  +  G  -  F)  x  n  =  0, 

(kxE  +  mXG-txF)  xn=0, 
since  according  to  (3-9)  they  are  transverse, 
k  •  E  =  m  •  G  =  ■{,  •  F=0. 

Using  the  experimental  coordinates  of  section  2,  one  has 

k  =  m  =  t  =0 
2  2  2 

and  a  system  of  6  complex  algebraic  equations  in  6  complex 
unknowns,  F  and  G. 


Before  solving  by  Cramer' s  rule  it  is  necessary  to  see  if  the 


system  determinant  has  any  significant  zeroes. 


1 

0 

0 

-1 

0 

0 1 

1 1 

-1 

0 

0 

0 

0 

0 

1 

0 

0 

-1 

0  : 

<3 

‘*3 

0 

0 

0 

0 

0 

^3 

0 

0 

-m3 

0 

0 

0 

1 

0 

-1 

0 

0  =  det 

^3 

0 

-“3 

0 

ml 

=  det  0 

0 

*3 

■^1 

-*3 

ml 

tl 

0 

*3 

0 

0 

0 

0 

0 

iQ 

0 

0 1 

io 

i 

0 

0 

mi 

0 

“3 

0 

0 

0 

0 

m3 
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Thus  the  system  is  block  diagonal  and  separates  into  two  inde¬ 
pendent  systems,  one  for  y  field  components  only  and  one  for 
x  and  z  components.  For  zeroes  of  the  y  system  determinant 


0  =  det 


-1 


-m 


=  l  -  m 

3  3 


or 


t  '  =  m  and  l  "  =  0 

3  3  3 


Using  (4.5)  the  latter  implies  l"  =  0  so  from  (4.10)  either 

oj  =  0  or  n  "  =  0.  In  addition 
o  1 


l  1  =  -V  cos^ 

3 

m  =  mcos7 

3 


so  from  (4.1),  (4.2),  (4.9),  (4.10)  said  cos0  =  -cos7 


n  1  cos$  =  cos0,  uj^O 
i  o 


=  n  1  Vl  -  sin20 


=  Vn  ' 2  -  sin2£ 
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L 


I 


which  implies  either  <u  =  0  or  n  ' 

*01 


=  1.  Neither  of  these  con¬ 
ditions,  m  =  0  or  n  '  =  1  and  n  "  =  0,  is  of  present  interest. 

r>  •»  i 


o  1  1 

For  the  x  and  z  system  determinant  to  vanish 


1 

c 

-1 

0 

0  =  det 

~ll 

■°3 

“l 

•tl 

1 3 

0 

0 

0 

0 

“l 

°3 

Implies  either  l 

3 

=  0 

=  m 

3 

or  m: 

=  t3m2  -  mgt2  (4.1?) 


former  is  of  no  interest.  Since  both  sides  of  the  latter  may 

be  varied  independently,  its  general  solution  is  that  both 

parts  are  equal  to  a  common  constant,  say  K.  For  the  left 

hand  side  then  K  is  real  and 

m2  =  m  2  +  m  2  =  Km 
13  3 

or 

n>12  +  Ua  -  *K)2  =  ¥K2, 


thus  possible  m^  are  |mj  s  ^K.  For  the  right  hand  side 


t2  =  t2  +  l2=Yll 
13  3 


and  l  is  complex  so 

3 


t  2  +  l 


I  2 


l  "2  =  Kt  ’  , 


(4.16  ) 
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at,  'i  "  =  Kt, " . 

3  3  3 

If  t  "  f  0  then  t  '  =  4k  and  t,  2  -  t,  1,2  =  rK2  so  possible  i 
3  3  13  1 

are  J-t,  |  2  From  (4.l6)  this  inequality  holds  even  if  l  " 

X  3 

The  only  solution  compatible  with  both  inequalities  and  (4.8) 
is 

\\\  -  &  -  UJ 

With  (4.16)  this  implies  t,  "  =  0,  the  no  absorption  case. 

3 

Then  from  (4.]),  (4.2),  (4.15),  and  (4.8) 


l  =  e  m 
3  13 

t,  =  m 


l  l 


and  eliminating  V  s  in  the  first  equation 


m2+e  2m2  =  e  (m  2  +  m  2) 

1  13  11  3 

or 


(e  -  l)m  2  +  (e  -  ®  2)m  2  -  0 

1  1  113 


26 


implying  e 


lorm^-eni^or 

1  13 


tan2e_  =  e  ,  e  " 

B  i  l 


=  0. 


(4-17) 


This  is  the  condition  for  Brewster's  angle.  For  polarization 
in  the  plane  of  incidence,  the  reflected  ray  vanishes  at  this 
angle.  It  has  been  shown  that  this  is  the  only  significant 
zero  of  the  system  determinant. 

Perhaps  this  investigation  of  null  solutions  is  only  a 
nicety  in  the  present  context.  It  does  endow  the  calculation 
with  a  satisfactory  feeling  of  completeness.  However  in  more 
complicated  cases  this  tool  will  still  apply  possibly  to  uncover 
regions  of  significant  departure  from  the  general  solution, 

--regions  which  may  be  as  useful  as  the  Brewster' s  angle  extinction. 

To  proceed  with  the  general  solution,  the  y  component 
system  gives  by  Cramer' s  rule 


k  -  ra 

F  =  -2 - 2.  e 


(4.18) 


m 


Likewise  for  the  x  component 


K  ra 

k  ra 

o  n 

F  =  -  E  , 

1  if  .  mf  1 

-t  m 


11 


(4.19) 


and  from  the  transversal ity  condition 


l 

F  ■  -  -si  F  .  (4.20) 

3  L  1 

3 

This  completes  the  determination  of  the  transmitted  fundamental . 

In  the  second  harmonic  boundary  value  problem  the  source 
wave  vector  is  given  by  (3.6)  and  is  complex. 

s  =  21  (4.21) 

I 

The  freely  propagating  wave  vector  magnitudes  are  from  (3-10) 

(4.22) 

(4.23) 

Continuity  of  fields  at  the  boundary  requires  equality  of 
phase  functions,  giving  for  the  complex  parts 

s"  •  x  =  t"  •  x,  n  •  x  =  0. 

But  by  (4.4)  and  (4.21)  _s"  has  no  tangential  part,  hence 

t  —  -  n* 


r  = 


4e  u)  2 
t2  =  -iL°- 
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As  before  all  wave  vectors  are  in  the  plane  of  incidence  and 


the  real  tangential  components  are  determined  by 


t' sinr  =  rsinp  =  s' sino. 


(' 


The  reflected  wave  vector  is  now  completely  determined  and  in 
particular  from  (4.21) 


cr  =  $ 


and  by  comparing  (4.24)  and  (4.8)  using  (4.22),  (4.21),  (4.2) 
and  (4.1) 


P 


=  0. 


0 


Now  define 


n  '  = 


ct' 

2tu  * 
o 


_  t.  Ct' 
n  =  ■— 
2  2d 


it 


(4.26) 


so  that 


=  n  • 2  -  n  "2  +  2tn  "  Vn  ' 2  -  sin20 
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Except  for  the  longitudinal  component  in  the  inhomogeneous 

term,  (4.28),  this  system  has  the  same  algebraic  form  as  before. 

For  the  y  component,  using  (4.29)  and  s  =  0, 

2 


(4.30) 


For  the  x  and  z  components 


r 

R  - - -  R  (4.31) 

a  r  x  \  i 

3 


where 


s  S 

3  1 


s  S 

1  3 


(s  P  -  s  P  ), 

3  1  13 


the  y  component  of  (4.29).  Solving  (4.28)  and  (4.29)  for  S^, 
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s 

x 


4tt  S  % 

—  (s  P  +  sP)  +  -4 - (s  P 

e  ii  33  Sse-e  31 
2  12 


s  P  ) 

1  3 


- - -  [(e  s2  -  e  s  2)p  .  c 

s2e  (e  -  e  )  2  11  1 

2  1  2 


and  using  (4.21),  (4.2),  (4.25),  and  (4.24) 


s  s  P  ] 

113  3 


S 

1 


4rre 

_ l 

s^e  (e  -  e  ) 

2  1  2 


-  S  8  P  ] 
13  3 


4n 


Lt  ^ 

3  1 


S  6  P  ]  . 

13  3 


Finally 


R 

l 


[t  P 

3  1 


+  t  P  ]. 

1  3 


(4.52) 
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5.  THE  REFLECTED  SECOND  HARMONIC  INTENSITY 


A  completely  polarized  beam  can  be  analyzed  by  finding 
the  intensitf^R  of  two  linearly  polarized  components  and  their 
phase  difference.  Arbitrarily  limiting  this  calculation  to  an 
examination  of  intensities  leaves  as  t.he  most  general  set  of 
measurements  the  four  polarization  cases  listed  in  section  2. 
Furthermore  all  four  cases  provide  different  information  from 
a  qualitative  standpoint.  Thus  in  the  last  section  the  inde¬ 
pendence  of  the  y  component  and  x  -  z  component  systems,  shows 
that  in  penetrating  the  surface,  a  ray  preserves  two  special 
cases  of  linear  polarization,  normal  or  parallel  to  the  plane 
of  incidence.  However  the  second  harmonic  polarization  process 
mixes  these  components  as  shown  by  (2.13)  and  (2.14). 


The  normal  component  second  harmonic  intensity  is  ,  from 


(4.30) 


with  P  generated  by  F  through  either  (2.13)  or  (2.14)  according 
2 

to  the  incident  polarization.  Thus 


2n  2  s  -  t  2 

J  =  9  - —  -3 - t3-  I v I  2  IF  I4  sin4Ycos2Y 

nn'e-€  r-t  '*>  1  2* 

12  3  3 


|x| 2  In2sin4Ycos2Y  (5-1) 
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where  the  second  J  index  denotes  incident  polarization.  Now 


I  =  IE  J2  =  |E  |2 
P  1  i‘  '3' 

so  by  the  tranversality  of  E 


Using  also  (2.14),  (4.19),  (4.1)  and  (4.2) 


|P  |2  =  £  |x|2  |F  2cosY(38in2¥  -  l)  +  F  2cosy|2 
2  1  3 

ni  k2  —  k  m2  JC  2 

=  t  lx!2  “2 - |G  I2  I  2cos2y|  1  +  -2-  ( 3 8in2Y  -2)|2, 

ml2  -  t  m2  1  P  l2 

3  3 

hence 


m  t2  -  l  m2 
3  3 


t  2 

x  Jl  +  -2_  (3sin2Y  -  2)  |  2cos2Y  (5-2) 

l 2 


Fran  (4.31)  and  (4.32) 
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JP  =■  lRJ"  -  Ir3!" 


2 

t  -  S 

3 

_ 3 _ 

e  -  e 

t  r2  -  r  t2 

1  2 

3 

3 

2 

It  P  +  t  P  I2 

1  3  1  13 


The  case  of  normal  component  illumination  simplifies  because  of 
the  absence  of  a  z  component  in  (2.13). 


J 

pn 


2rr 

2 

t  -  S 

3  3 

e 

-  e 

t  r2  -  r  t2 

l 

2 

3  3 

2n 

2 

t  -  6 

3 

3  1 

e  -  e 

t  r2  -  r  t2 

1  2 

3 

3  1 

2 


2 


It  xF  £|2sin2Y(3cos2Y  -  l)2 

S  2 


2 


Ixl2  * 


x  In2sin2Y(3cos2Y  -  l)2 


(5.5) 


For  parallel  component  illumination,  from  (2.14) 


2 


P  =  \  x( -5F  2sin¥cos2¥  +  F  2sinY) 

11  3 


-i 


/  m  k2  -  k  m: 


2\  2 


lit2  •  In21 


[1  -  -2-  (3cossY  < 


P  =vFF  sinY 

3  13 


/m  k2  -  k  m2\  It 
3  3  1  E  2  — i-3.  sinY 

m  l2  -  l  ms  I  1  k2 

3  3  1  3 


hence 


JPP 


at® 

2 

t  -  S 

m  k2  -  k  m2 

3  3 

i 

3  3 

e  -  c 

t  r2  -  r  t2 

1  t2  ■  t  B2 

1  2 

'3  3 

3  3 

Ixl 


t  [1 - (3cos2¥  +  l)] 


1 1 
2t  -ua. 

1  t2 


l) ] sinY 


2I  2sin2Y 
P 


(5.4) 
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6.  NUMERICAL  RESULTS 


The  above  intensities  are  functions  of  a  mmber  of  parameters, 

J  =  JJk'<V  *1*  e2*  \*  9> 

Also  several  factors  are  complex  algebraic  so  one  may  feel  it  is 
hard  to  extract  information  at  a  glance.  For  discussion  purposes 
define A4 

t  2  T  i  l 

f( 0,  Y)  =  1 1  — 3—  (3cos2Y  +  1)  -  2  -i  |  2sin2Y  «  j  (y)  (6.1) 

l2  T  t2  55 

3 


l  2 

g(0,  Y)  =  |l  +  —  (3ain2Y  -  2) |  2cos2Y  <*  J  (y)  (6.2) 

These  two  factors  contain  all  the  Y  dependence  of  two  of  the  inten¬ 
sities  and 

f(0,  \  -  Y)  <*  JjY)  (6.3) 

g(0,  \  -  Y)  <x  jpn(Y)  (6.4) 

contain  the  Y  dependence  of  the  other  two.  By  changing  the 
definition  of  Y  from  a  geometrical  one  to  a  physical  one  it  is 
possible  to  ascribe  major  features  of  the  variation  of  intensity 

14  The  symbol  «  indicates  proportionality. 
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with  crystal  orientation  to  crystalline  features.  Let  V  be  the 
angle  from  [100]  to  the  projection  of  the  incident  electric  vector 
on  the  reflecting  surface.  Then  in  (6.5)  and  (6.4)  ¥’  =  |  -  Y 
since  the  incident  polarization  is  normal  but  in  (6.1)  and  (6.2) 

V  =  Y.  The  curve  f(0,  V  )  shows  extinctions  on  the  [100]  and 
[110]  axes  and  a  maximun  on  the  [ill]  axis  (V  =  cos"^"^  =  54.7*). 
The  curve  g(0,  ¥’ )  shows  extinctions  at  the  [100]  and  [ill]  axes 
and  an  absolute  maxinnxn  at  the  [110]  axis.  These  two  curves 
(see  Pig.  5)  are  independent  of  0  and  the  material  constants. 
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Using  the  frequency  and  material  constants  given  in  section  2, 
f(e,  Y)  and  g( 9,  Y)  have  been  computed  for  a  representative 
sampling  of  points  (0,  Y) .  The  comparison  of  f(0,  Y)  and  f(0,  Y) 
may  be  regarded  as  exhibiting  the  difference  in  shape  between 
J^Y’ )  an<*  Jpp^*  Likewise  comparison  of  g(0,  Y)  to  g(0,  Y) 
exhibits  tne  difference  between  J  (Y’)  and  J  (Y).  Since  the  Y 
dependence  is  arbitrary  to  within  a  constant  multiple,  the  com¬ 
parison  in  shape  is  made  by  normalizing  the  maxima.  For  8  =  45 

the  maximun  variation  in  shape  between  J  (Y')  and  J  (Y)  is 

on  pp 

60j6  based  on  f(0,  Y)  and  the  maximum  variation  in  shape  between 
J  (Y')  and  J  (Y)  is  14$.  tn  both  cases  the  maximum  variation 
occurs  when  the  projection  of  the  incident  electric  vector  falls 
between  the  [100]  and  [Hi]  axes.  For  larger  0  the  differences 
increase.  Thus  though  J^(  Y)  and  J  ( y)  look  more  complicated,  they  have 
much  the  same  shape  as  J^Y' )  and  J  (  Y' )  respectively.  As  to 
specific  details  of  difference  J^(  y)  does  not  show  a  perfect 
extinction  at  [110]  and  neither  does  Y)  at  [ill].  Furthermore 
the  maximum  of  J  (Y)  and  the  relative  minimum  of  J  (y)  show  a 
small  shift  from  [111] . 

In  addition  f  and  g  were  computed  in  the  approximations 

h  —  0  and  h  -  0  s  k  which  are  convenient  since  they  allow 
112 

partial  or  complete  shift  to  real  arithmetic.  At  first  glance 

the  last  approximation  seems  poor  since  it  means  neglecting 

h  =  1.53  by  comparison  with  v  =  3.89.  The  first  approximation 
2  2 

changes  f(45*,  Y)  by  .05$  near  its  maximum  and  the  second  by 
. 01$  at  the  same  point.  In  these  approximations  the  largest 
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changes  relative  to  the  exact  values  at  the  same  points  are  10.3$ 

at  Y  =  90*  and  24$  near  ¥  =  85*.  Both  approximations  are  the  same 

for  g(45*,  ¥)  and  change  it  by  .2$.  It  seems  clear  that  once  a 

very  small  absorption  is  admitted,  increasing  the  absorption 

causes  no  qualitative  change  until  the  critical  point  n"  =  n’  is 

reached,  at  which  point  propagation  ceases.  The  latter  is  vaguely 

indicated  above  by  failure  of  the  analysis.  Equation  (4.13) 

requires  n  '  >  n  "  unless  e  '  £  0  and  the  latter  causes  a  funda- 
111 

mental  change  in  the  character  of  the  Helmhotz  equation  (3 .4). 
Furthermore,  even  the  transition  from  no  absorption  to  some  absorp¬ 
tion  appears  smooth.  The  only  qualitative  change  indicated  above 
is  the  disappearance  of  the  system  null  solution  (4.17).  Evidently 
when  there  is  seme  absorption,  the  Brewster  extinction  becomes 
only  a  relative  minimum. 

The  next  most  complicated  factors  in  the  intensities  are 
recognizable  as  transmission  ratios — one  for  transmission  of  the 
fundamental  into  the  medium  and  one  for  transmission  of  the  second 
harmonic  out.  The  exact  evaluation  of  the  intensity  coefficients 


for  0  =  45*  gives  the  following. 

=  2.926  X  10’4|xl2In2f(0,  ¥•)  (6.5) 

Jnp  =  6.279  X  10"4|x|2ip2g(45*,  V)  (6.6) 

Jpn  =  4,548  x  10’4lxi2In26(0>  *’)  (6.7) 

Jpp  =  9.529  x  10’4|x|  2Ip2T(k5*>  V)  (6.8) 


40 


The  following  table  lists  per  cent  errors  in  the  numerical 
coefficients  in  the  approximations  indicated. 


J  J  J  J 

nn  np  pn  pp 


Table  2,  Errors  in  Two  Approximations  of  the  Intensities 

The  conclusions  of  the  last  paragraph  on  the  general  effects  of 

absorption  still  apply.  Thus  absorption  considerations  can  be 

neglected  if  only  qualitative  behavior  is  sought. 

Much  of  the  behavior  expected  has  been  verified  by  experiment. 

The  variation  of  J  and  J  with  Y  on  an  arbitrary  intensity 
nn  pn 

scale  were  reported  in  1961  by  Ducuing  and  Bloembergen. 15  Certain 

aspects  of  the  variation  in  $  have  been  obtained  including  the 

minimum  near  the  Brewster  extinction. 16  Finally  |x|  relative  to 

that  of  potassium  dihydrogen  phosphate  has  been  obtained  at  several 

frequencies  and  for  other  semiconductors  and  the  phase  of  x  for 

GaAs  has  been  obtained  at  ruby  laser  frequency.17  The  data  on 

|x|  is  sufficient  to  show  its  general  dispersion  characteristics 

NL  L 

and  it  is  seen  that  peaks  of  x  (2tu)  roughly  match  those  in  x  (tu) . 

15  See  footnote  5.  Hie  same  results  appear  in  N.  Bloembergen, 
Nonlinear  Optics,  Benjamin  1965,  PP*  126-7  with  better  identification 
of  polarization  cases. 

16  R.  K.  Chang  and  N.  Bloembergen,  Phys.  Rev.  144,  775  (1966) 

17  R.  K.  Chang,  J.  Ducuing  ana  N.  Bloembergen,  Phys.  Rev.  Letters  15, 
415  (1965).  For  a  review  of  reflected  second  harmonic  experiments, 
see  N.  Bloembergen,  Optica  Acta  1J,  311  (1966). 


7 .  REMARKS  ON  THE  GENERAL  PHASE  MATCHING  PROBLEM 


In  the  first  section  the  possibility  that  the  output  signal 
might  attain  intensities  of  the  same  order  of  magnitude  as  the 
pump  signal  was  mentioned.  When  the  output  signal  is  the  second 
harmonic,  phase  match  is  known  to  be  a  necessary  condition  for 
this  possibility.18  In  this  exposition  of  second  harmonic  genera¬ 
tion  by  reflection,  phase  match  occurred  as  a  null  solution  to 

the  nonlinear  wave  equation,  i.e.,  e  =  e  in  (^.8).  Because  the 

l  2 

concept  of  null  solution  is  quite  general,  there  is  a  strong 
inference  that  it  is  an  avenue  for  extending  those  implications 
of  phase  matching  which  are  contained  in  the  above  necessary 
condition.  Thus  it  should  be  possible  to  generalize  the  results 
listed  on  pp.  1-2  for  anisotropic  nonabsorbing  media  and  for 
optically  active  nonab sorbing  media  to  anisotropic  media  exhibiting 
the  most  general  linear  optical  effects. 

The  character  of  this  generalization  may  be  sketched  in  outline. 
It  will  be  matrix  algebraic  and  so  will  not  differ  much  from 
Kleinman’ s  approach.1  If  the  indications  of  the  last  section  on 
the  relative  unimportance  of  absorptive  effects  hold,  the  generaliza¬ 
tion  of  the  work  on  birefringent  media  to  include  anisotropic 
absorption  or  the  generalization  of  the  work  on  optical  activity 
to  include  circular  dichroism  is  trivial.  However,  as  a  survey 
of  effects  attributable  to  the  imaginary  part  of  the  dielectric 

18  J.  A.  Armstrong,  N,  Bloembergen,  J.  Ducuing  and  P.  S.  Pershan 
Phys.  Rev.  1?7,  19^0  (1962). 
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constant,  the  present  effort  is  incomplete  due  partly  to  the  ex¬ 
clusion  of  total  reflection  and  partly  to  the  lack  of  special 
consideration  of  the  relation  between  optical  activity  \nd  circular 
dichroism.  Thus  it  will  be  necessary  to  complete  the  description 
of  the  dielectric  constant.  These  generalizations  still  leave 
two  classes  of  optical  effects  disjoint  under  the  capability  of 
describing  phase  match  in  terms  of  material  properties.  The  two 
classes  could  be  called  symmetric  and  antisymmetric.19  Since  the 
dielectric  constant  for  optical  activity  requires  wave  vector 
dependence,  perhaps  these  two  classes  can  be  merged  by  generalizing 
the  anzatz  (2.4)  to  include  spatial  dispersion. 
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19  See  the  categorization  of  linear  optics,  p.  2. 
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